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We study a uniform matrix product state as a variational state for classical and quantum spin 
chains in the thermodynamic limit. Under a careful treatment of the translational symmetry, 
eigen values of the transfer matrix defined in the calculation of expectation values can reflect 
the periodicity of the ground state and indicate optimum periodicity of the matrix product 
state. We discuss the relation between the periodicity and accuracy of magnetization curves. 
This approach is free from the error due to finite system size, which works well especially for 
the magnetic plateau problem. 
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A recent trend of numerical calculations of quantum 
mechanics is a construction of new variational states as 
extensions of the matrix product state (MPS), such as 
the tensor product state (TPS),^'^-' the projected entan- 
gled pair state (PEPS),"^^ the tree tensor network state 
(TTNS),^'^-* and the multi-scale entanglement renormal- 
ization ansatz (MERA) state. The most driving force of 
this trend is the success of the density matrix renormal- 
ization group (DMRG) method, which is a power- 
ful numerical method for one-dimensional strongly cor- 
related systems. The MPS is not only used as a vari- 
ational state^^^ but also the exact ground state of the 
valence-bond-solid (VBS) state in the Affleck-Kennedy- 
Lieb-Tasaki (AKLT) model, and the exact eigen states 
of the Bethe ansatz.". is) 

The extensions of the MPS in the thermodynamic 
limit have been widely studied, where we expect that 
the MPS has uniform, i.e., site independent, matrices 
for the uniform Hamiltonian, for example, in the infi- 
nite time-evolving block decimation (iTEBD), the in- 
finite PEPS (iPEPS),") and scale invariant MERA.^^) 
Such a uniform matrix is obtained as a fixed point of the 
DMRG.i°'") As a merit of the uniform MPS, we can 
handle the infinite uniform system directly, that is, we 
can construct variational states with finite dimensional 
matrices in spite of the infinite dimensional Hilbert space 
of the infinite system. In some methods, the matrices are 
optimized by the iterative projection with using Suzuki- 
Trotter checker-board decomposition.^^' 

A natural question for the uniform MPS is how to deal 
with the spontaneous symmetry breaking of the transla- 
tional symmetry, where a ground state has long peri- 
odicity. One of interesting examples is the magnetiza- 
tion plateau, which emerges in some frustrated systems 
and has been widely researched, because we can expect 
a new quantum phase in a new plateau. For example, in 
the 5 = 1/2 antiferromagnetic zigzag chain, the novel 
even-odd effect of the 1/3 plateau has been studied by 
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using the DMRG.^^' In the calculation of the magnetiza- 
tion curve, the finite size effect causes the small steps in 
the magnetization. On the contrary, the uniform MPS in 
the thermodynamic limit gives a smooth magnetization 
curve by definition, which is an advantage of the uniform 
MPS. 

In this letter, we study the uniform MPS with arbi- 
trary periodicity and a boundary matrix for magnetic 
plateaus. To treat the translational symmetry explicitly, 
we do not use the checker-board decomposition, which 
breaks the primitive periodicity, i.e., the one-site trans- 
lational symmetry. When we consider a linear combina- 
tion of g-fold degenerated ground states with (/-site peri- 
odicity, the periodicity of the MPS is determined by the 
boundary matrix. The boundary matrix is important also 
in discussions of the fixed point of the DMRG^°' ^^-^ and in 
the continuous MPS.^^^ The MPS used in this letter has 
matrices generalized toward Baxter's interaction-round- 
a-face(IRF)-type MPS,^'^-' which has an advantage in the 
case of dimerization. Using the IRF-type MPS as a gen- 
eralization of the usual MPS, we demonstrate the perfor- 
mance to calculate the magnetization curve in quantum 
spin chains. From the numerical point of view, it is im- 
portant to consider the optimum periodicity of the MPS 
in order to obtain the best variational ground state. 

Let us start from the definition of a finite iV-site vari- 
ational state 1 4') with the boundary matrix Aq. The 
total Hilbert space for a spin S system is spanned by 
|cr) = rii^r I'''*) '^i*^ (2*^ + l)-dimensional local basis 
\<Ji), where cr — {ai, . . . , ctat}. Many theoretical stud- 
ies have tried to reduce computational memory by us- 
ing tensors A with artificial local bases a,; instead of 
(25" -I- 1)^ dimensional vector \l/(cr) = {cr\^). Gener- 
ally, the number of tensors A and local bases ai de- 
pend on the coordination number and the network of 
tensors. In one dimensional case, a simple MPS with- 
out Ao is defined as *(cr) = Ea Il^i ^?a,,a,+i . where 
a = {ai, . . . , aAf}, and a^+i — ai- We suppose the di- 
mension of each is finite m. At m = {2S -\- 1)^ any 
^ will be expressed. There are other expressions such 
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Fig. 1. Graphical representations of the usual MPS and the IRF- 
type MPS. Filling the small square means contraction with re- 
spect to the local artificial bases Qi. 
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l^-) = Tr 
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with {A,)a.,a' = J2a 

As a generalization of the usual MPS, we use a IRF- 
type MPS 



(Tjv.cri 
N ,a 



(1) 



where a — {ai,...,ajv} are artificial local bases, and 
Aq = An is a generalized boundary matrix. We call 
it a IRF-type MPS in distinction from the usual MPS 
with the boundary matrix,^"' which is represented 
as the special case that A^^^, = Af.^^, and A^^^^, = 
Sa" ^Ar Qa"^a"a'- The IRF-typc MPS is an extension 
of the usual MPS as shown in Fig. 1, and can represent 
local entanglement very well as shown below. 

To demonstrate the ability of the MPS at small m 
we consider the direct product state of the local singlets 
which is the exact ground state of the spin-1/2 dimer- 
ized Hamiltonian Hd = Si ^zi-i • 82^. Corresponding 
the two-site translational invariance of Hu , we suppose 
A21-1 = Ai, and ^2* = ^2 for all i. The usual MPS at 
m = 1 is just a direct product state of the local spins, 
which is a classical state with no quantum entanglement. 
At least, m = 2 is required to represent the exact ground 
state for the usual MPS.^^) The IRF-type MPS can rep- 
resent the exact ground state even in the m = 1 case.^^^ 
Here, Aq is identical and is not important. 

To show the importance of Aq, let us consider the 
doubly degenerated classical Neel states expressed by 
the IRF-type MPS, which is realized at m = 1 by 



j^n = ^ and AJ^ = Aj^ = 0. Then, we ob- 
tain l^-) A^l^^^l n ■■■) + it • • • > for even N. 
We emphasize the boundary matrix Aq"^ determines the 
periodicity of j^*), and l^*) has one-site translational in- 
variance if A^'^ = AI^. 

In general, there is the exact mapping from the IRF- 
type MPS with m = tuq to the usual MPS with 
m = (25 + l)mo via the singular value decomposi- 
tion (SVD). If we apply the SVD to IRF-type MPS 
as Af-;^, = J:p{'^c,\U^\^HpmV^H^'a') with singu- 



lar values 

mension of the index /? is the 



we obtain the usual MPS A, 



where the di- 
as era, that is. 



same as era 
{2S + l)mQ. 

Let us summarize the formulations of the IRF-type 
MPS with the finite period p and the dimension m to 
deal with the thermodynamic limit N ^ 00. We suppose 
translational invariance with period p as Ai = Ai^np for 
all non-negative integer i and n. The system size has the 
relation N — Nip + Nq with the non-negative integer 
A^i and Nq — N mod p. We consider the mismatch A'o 
generally, but the effect due to Nq will be merged into the 
boundary effect due to Aq as shown later. For simplicity, 
we consider Aq = Ap and Nq = for a while. 

For a finite N, the inner product is written by (^'l^') — 



Tr 



rN 



n.=i 

ments 



T 



with a transfer matrix T,- with matrix ele- 



i\OL" ol'" ■ 



(2) 



Using the p times periodicity for the N — N\p system, 
we have (^-l^-) = Tr [T^^] with T = OLi ■ 

In the thermodynamic limit, only principal eigenval- 
ues Xj of T is important, where "principal" means that 
|Aj| = A is maximum. Hereafter, we suppose T is diago- 
nalizable. We denote principal eigenvalues as Xj = X e*^'' 
and corresponding left (and right) eigenvector as Uj (vj) 



with the bi-orthogonal condition v]uj 



Ni, we have 



At large 



(3) 



Hereafter, we renormalize Ai = X'^p Ai and rewrite it as 
Ai. Then, we obtain lim7v-i.oo(^|^) — "^j^j with nor- 
malization Vj — limjv-i.00 e^^J^i. By definition, N ^ 00 
means A^i — > 00. For a general phase 9, we can adopt 



lim 



iBNi 



,9^0 
,9^0 



(4) 



However, a normalization i^j for a commensurate phase 
9j — 2-Kj/q is able to be non-zero if we consider the 
special case Ni ex: q. This means that emergence of the 
commensurate phase indicates gp-periodicity of the state. 

To control the periodicity in the thermodynamic limit, 
we use the boundary operator ^o- We can define the 
boundary transfer matrix Tq for p + Nq sites including 
Aq and the mismatch Nq. Then, the norm is written as 
(5f|5-) = Tr [TqT^^-^] = J^j with 



lim 



N- 



(5) 



for TV — Nip + Nq. Here, whether Vj is zero or not is 
finally determined by the choice of the boundary transfer 
matrix Tq. Since one of principal eigenvalues of T is real, 
which will be shown in our numerical calculations, we 
will choose Tq to eliminate contribution from the non- 
real principal eigenvalue and to set Vj = 5ji supposing 
9i = 0, which means the p-site translational invariance 
and enables us to express the expectation values in a 
simple form. 

For the calculation of the bulk energy e — 
limAT^oo jvl^^ l ^) ' '^an neglect the boundary effect 
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Fig. 2. Magnetization per site M and principal eigenvalues of the 
matrix T as a function of the model parameter t in the Hamil- 
tonian 'Hj(t). 

which gives the 0{1/N) boundary energy. Here the local 
Hamiltonian Hii^i'] are supposed to include i' — i + 1 spins 
from ith site to z'th site. Then, after using Vj = Sji^ we 
obtain e — v|/iui with 

naa"CTi;a'a"Vjv^p - Nhp - 2 

Where IM^^^ = k> (nfjT' ^^""0 . 

\ / act' 

(T = (72, . . . ,a-NhP-i- Here iV;i is supposed to be enough 
large to include the periodicity of the Hamiltonian. In 
the same manner, one can calculate any expectation val- 
ues per site. To optimize the IRF-type MPS the down- 
hill simplex method and the modified Powell method are 
used.26) 

To reveal the relation between principal eigenvalues Aj 
with p — 1 and the spontaneous symmetry breaking of 
the translational symmetry, first, we discuss the magneti- 
zation plateau state of the classical S — 1/2 Ising model 
with the third next nearest neighbor interaction given by 

ni{t) = {mS! + ELi Jkit)S!S^^^) , where Sf 
represents the z component of the spin-1/2 operator. The 
longitudinal magnetic field and fc-th nearest neighbor in- 
teraction are given by H and J| respectively. To show 
various magnetic phases of the ground state, all param- 
eters are defined as a function of one parameter t E [0,1] 
as follows; H{t) = 0.9i, Jf(i) = -0.14t 0.57, J|(i) = 
-0.22i -I- 0.46, J§{t) = -0.48t + 0.38. In fact, we detect 
six different phases in the parameter region of t with us- 
ing analytic results. 

Figure 2 shows the magnetization per site as a func- 
tion of t obtained by the IRF-type MPS with p — 1 and 
m — 3, which completely agrees with that of the ana- 
lytic results. The spin configuration from analytic results 
in each phase are also denoted by the black allow. The 
six diagrams mean the complex principal eigenvalues Xj 
of the matrix T in each phase. It is shown that q princi- 
pal eigenvalues correspond to the spontaneous symmetry 
breaking of the g-site translation. If we enlarge the peri- 
odicity p to g, analytic results can be exactly represented 
at TO = 1 . In this case, the unique principal eigenvalue is 
obtained. 
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Fig. 3. Magnetization per site M as a function of magnetic field 
H for (a)S' = 1/2 and (b)S = 1 Heisenberg chains with the 
periodicity p and dimension m. 

Secondly, we show the result of the uniform quan- 
tum chain. The Hamiltonian of the Heisenberg chain 
with the longitudinal magnetic field is given by Hhb = 
J2i {HS^ + JSj • Sj_,_;i), where S^ represents the spin-S' 
operator. Here, we take S* = 1/2 and 1. The parameter 
J means the intensity of the isotropic exchange inter- 
action, and it is taken as the unit of energy. The MPS 
with the period p = 2 is selected in both spin-S" chains, 
because the MPS with p — 2 gives the minimum energy 
for fixed m within p < 6, which reminds us Neel state's 
periodicity. 

Figure 3 shows the magnetization curve for (a) S — 
1/2 and (b) S = 1 Heisenberg chains. The dimension of 
the matrix to is up to three in each chain. As reference 
data, we also show (a) the rigorous result in S" = 1/2 
from the Bethe ansatz, and (b) the reported data from 
the DMRG calculation, respectively. The results of our 
MPS of TO = 3 agree with the references in this scale. 

Finally, we show the result oi S = 1/2 Ising-like 
zigzag XXZ chain, which is a uniform quantum frus- 
trated chain. The Hamiltonian is given by "Hxxz = 
E,(if5f+ELi-/|(A(5f5f+,+SfSf+,)+SfSf+,)),with 
the anisotropic parameter A. We choose parameters as 
A = 0.1, Jf = 0.4, and J| = 1.0. It is known that the 
1/3-plateau state appears in this system under the pa- 
rameters. In addition, q — i principal eigenvalues are 
confirmed in the matrix T with the period p = 1 in the 
1/3-plateau. The DMRG data shown by the black line 
in Fig. 4 are also re-calculated following Ref. 21). The 
perfect 1/3 plateau is obtained by the MPS in small m, 
while there is the finite-size effect for the DMRG result. 
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Fig. 4. Magnetization per site M as a function of magnetic field 
H for the S = 1/2 zigzag chain "Hxxz- The DMRG data is 
calculated for 192 sites with the open boundary condition. The 
dotted line means the magnetization M of the best variational 
state. 



In the plateau we obtain the minimum energy at p = 3. 
Unhke Hhb, "Hxxz shows the transition from p — i io 
p = 4 in large magnetic fields as shown in Fig. 4. The 
transition is artificial due to small m, but the DMRG 
data also show the property change about the stability 
of the odd M states. It can be concluded that more large 
m is needed. 

In summary, we investigated the uniform IRF-type 
MPS with an arbitrary period p and with a boundary ma- 
trix as a generalization of the usual MPS. The boundary 
matrix related to the periodicity is used to hold the p- 
site translational invariance in the numerical calculation 
for the classical or quantum spin chain in the thermody- 
namic limit. We found that the number of the principal 
eigenvalues, 9, in the matrix T with p 1 is consistent 
with the period of the spontaneous symmetry braking 
ground state in the classical Ising model and indicates the 
optimum periodicity which gives the minimum energy for 
the fixed dimension of the artificial basis a. The IRF-type 
MPS with the optimum periodicity p = 2 gives consistent 
magnetization curves in 5* = 1/2 and 1 Heisenberg chains 
and that with p = 2> captures the correct magnetization 
at the 1/3 plateau without the finite-size effect in the 
S* = 1/2 zigzag chain. This means that the plateau state 
is easy to reproduce via this variational state which has 
small quantum entanglement except for around the up- 
per critical magnetic field of the 1/3 plateau. Our demon- 
stration is limited for small m, but smallness of m leads 
us to the optimum periodicity p = 4 shown in the mag- 
netization curve in zigzag chain, which can be related to 
property of the true ground state. To establish an opti- 
mization method for large m is one of future works. 

Note that the period of the MPS is important for the 
wave function prediction in the infinite algorithm of the 
DMRG method efficiently.^®"^^^ The number of principal 
eigenvalues in the matrix T is also meaningful informa- 
tion for the prediction. 
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